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Abstract 

We solve Schrodinger's equation for the ground-state of /owr-dimensional Yang-Mills theory 
as an expansion in inverse powers of the coupling. Expectation values computed with the 
leading order approximation are reduced to a calculation in two-dimensional Yang-Mills 
theory which is known to confine. Consequently the Wilson loop in the four-dimensional 
theory obeys an area law to leading order and the coupling becomes infinite as the mass- 
scale goes to zero. 
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Continuum Strong Coupling Expansion of Yang-Mills Theory: 
Quark Confinement and Infra-Red Slavery 



1. Introduction 



Despite considerable success in describing high energy phenomena it remains to be 
shown that Quantum Chromo dynamics accounts for the known and conjectured properties 
of the Strong Interactions at low energies such as quark confinement and chiral symmetry 
breakdown. The strong coupling expansion on the lattice [1], which leads to an area law 
for Wilson loops, provides a framework for the understanding of confinement but only 
in a calculational regime where the continuous nature of space-time has been spoilt and 
the connection with continuum field theory has only been established numerically. The 
difficulty with attempting a strong-coupling expansion directly in the continuum is that 
our conventional understanding of quantum field theory is based on a different kind of 
expansion, namely the semi-classical, or loop, expansion. Indeed, the particle interpreta- 
tion of field theory derives from analysing the free part of the action in terms of harmonic 
oscillators. This is encoded into the LSZ formalism which underpins the use of field theory 
in Particle Physics. In Yang-Mills theory asymptotic freedom [2] implies that the LSZ 
realisation of field theory is applicable at high energies. As the energy scale is reduced the 
running coupling grows but it is not known if it continues to grow beyond the perturbative 
regime, although it is usually thought to do so, giving rise to infra-red slavery. In ordinary 
Quantum Mechanics the Born approximation, in contrast to the semi-classical approxima- 
tion, provides an expansion in decreasing powers of ft 2 . In this paper we develop such a 
scheme for continuum Yang-Mills theory in four dimensions. We will solve Schrodinger's 
equation for the wave-functional of the ground-state as an expansion in inverse powers 
of the coupling. Our strategy will be based on an ansatz in which we assume that the 
logarithm of the wave-functional is the integral of a local function of the gauge-potential. 
This enables us to relate ground-state expectation values in the four-dimensional theory 
to calculations in two-dimensional Yang-Mills theory. It is well known that Yang-Mills 
theory in two dimensions confines, and as a consequence, our leading order approximation 
for the Wilson loop in the four-dimensional theory has area-law behaviour. Furthermore 
the coupling depends on an arbitrary mass-scale, //, growing stronger as \i is reduced. 
The relationship between two-dimensional Yang-Mills theory and certain string theories 
has received attention recently [3], so our result provides a link between this and the 
four-dimensional theory. 

The Yang-Mills action can be written with the coupling factored out 



(The gauge potentials are elements of the Lie algebra of the gauge group so A = A C T C \ 



and (T C Y = -T c , tr (T A T B ) = -5 AB , [T A ,T B ] = fABC T c^ g 2 playg the role of 



Planck's constant since the action occurs in the path-integral as expiS/h . Thus the LSZ 




E = -A + VA + [A, A ], B = VAA + AAA. 



(1.1) 
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realisation of Yang-Mills theory is appropriate to the semi-classical regime where physical 
quantities are expandable as power series in g 2 . If we want to understand the strong 
coupling limit of the theory we must expand in powers of 1 / g 2 , in other words we need to 
study the ultra-quantum mechanical limit where the (essentially semi-classical) notion of 
particles may not be relevant. In this paper we will solve Schrodinger's equation using such 
an expansion for the ground-state. Of course, either kind of expansion should be viewed as 
being of a formal nature since renormalisation leads to the introduction of a length scale 
in terms of which the running coupling is fixed as a function of a mass scale. Rather, the 
two formal power series should be viewed as different ways of organising the calculation. 

To simplify things let's consider Schrodinger's equation for a quantum mechanical 
particle of mass m moving in a potential V(x) 

h 2 d 2 \ 

+ V(x) # = (1.2) 



2m dx 2 

If we write the wave-function as \1/ = exp — W(x) then this becomes 



h 2 



2m 



((W') 2 -W")+V = £. (1.3) 



In the semi-classical approximation, which corresponds to the loop expansion in the usual 
perturbative approach to field theory, we expand 

oo 

W = J2h n - 1 W n , (1.4) 
so that to leading order we can neglect W" and obtain the Hamilton-Jacobi equation 

-^Wf + v = £, <"> 

which identifies Wo as the classical action. When the potential is sufficently weak we can 
use instead the Born approximation. Assuming that V — > as |x| — > oo and that there are 
no bound states implies that the lowest value of S is zero corresponding to a (non-square- 
integrable) wave-function that is constant at spatial infinity. If we take \I/ = 1 — ip, with ip 
a small perturbation, then re-arranging (1.2) gives 

id 2 y\ _1 y 



^={-^ + ¥) ¥ (L6) 

which can be expanded in powers of V/Ti 2 . 

Another way of organising the same expansion is to write W as a power series in 1/h 2 , 
W = Y^=i n W n and expand (1.3). To leading order and next to leading order this 
gives (with S = 0) 

J-W» + V = 0, W£ = (W{) 2 . (1.7) 
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It is easy to check that this gives the same results to this order as expanding (1.6). This 
is the approach we will take to the Schrodinger equation for Yang-Mills theory. 

We will follow the formulation of the Schrodinger equation for Yang-Mills given by 
Jackiw in [4], which has the merits of being both simple and closely related to ordinary 
quantum mechanics. (This reference also explains the modification to the wave-functional 
required when there is a non-zero 6>-angle.) We work in the Weyl gauge (Aq = 0) so that 
the canonical coordinates are A c (x) and their conjugate momenta are — g~ 2 ~E c (x.). The 
Hamiltonian is 

#[A,E] = J d 3 xtr (E 2 + B 2 ) . (1.8) 

Since A (x) does not appear in the action the Euler-Lagrange equation obtained by varying 
A)(x), 

V-E + A-E-E-A = 0, (1.9) 

is a constraint rather than a Hamiltonian equation of motion. In the Schrodinger repre- 
sentation with % = 1 E c is represented by the operator 

EC < X > = ^nfej' <L10) 
(showing explicitly that g 2 plays the role of Planck's constant). Schrodinger's equation is 

-i 9 2 A + 9 - 2 B)*[A] = f<P[A], 
and the constraint is imposed as 

Now T is the generator of infinitesimal gauge transformations acting on functionals of A so 
this constraint simply implies that the wave- functional is gauge invariant. The Schrodinger 
equation stands in need of regularisation because the two functional derivatives act at the 
same point in space. In the next section we shall construct a regulated operator, A^, 
that preserves the gauge invariance and underlying geometry of the Hamiltonian operator. 
It acts on integrals of local functions of A to produce other local integrals with finite 
coefficents, unlike A which produces coefficents that are ill-defined. This means that we 
have a Hamiltonian operator that has a finite action on the space of sums of products of 
integrals of local functions of A with finite coefficents. The price for being able to do this 
is, as always, the introduction of an arbitrariness in the renormalisation procedure that, 
in this case, is encoded into a function, 0. This is needed to separate the finite part from 
quantities that diverge as we remove the cut-off. If the theory is renormalisable, as we 



4 



hope, we will be able to absorb this arbitrariness into the coupling constant so that physical 
quantities will not depend on 0. Thus the finite version of the Schrodinger equation will 
be 

{-\g 2 ^ + 4^) *[A] = £*[A] (1.13) 

where g is now a functional of 4>. If we attempt a strong-coupling expansion of the ground- 
state wave-functional 



tt[A] = exp - S[A] =expl-J2 9^ n S n [A] , S = ^ (1.14) 

\ n=l / n=l 

the Schrodinger equation becomes 

I^ S -^(/^._^) + _l B = £ (1 , 5) 



so that to lowest order 



A 4> S 1 = 2{£ 1 -B). (1.16) 



The right-hand side is the integral of a gauge invariant function that is local in the sense 
of being constructed from fields at the point x and a finite number of derivatives. By 
construction, A<f, will act on integrals of local gauge-invariant functions to produce other 
integrals of local gauge-invariant fuctions with finite coefficents. So we can look for a 
solution for S\ that is itself the integral of a local function. In the next section we will 
show that 

A^B = Pi( [ d 3 x)+[3 2 B, (1.17) 



where, for the gauge group SU(N) the coefficents are ft = (p^(0) (N 2 - l)/32O0r 3 , and 
ft = -110' (0) N/12^ 3 . So if we take 

Si = ~ £i = -w f d " x ^ (i-is) 

P2 P2 J 

then equation (1.16) is satisfied. The infra-red divergent integral over space in S\ is what 
we would expect for the energy of the ground-state, since by translational invariance the 
energy density must be constant. More generally we adopt an ansatz of locality as a 
strategy for solving the Schrodinger equation in a strong coupling expansion. We will 
assume that S[A] is a sum of integrals of local functions of A. Then, by construction 
A^S will also be a local integral, and so too will the other terms entering Schrodinger's 
equation in the form (1.15). Taking S to be gauge-invariant satisfies the constraint (1.12), 
and taking the integrals to be invariant under spatial translations and rotations acting on 
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A ensures that we are studying the ground-state which can carry neither momentum nor 
angular momentum. We also require that the ground-state wave-functional be Lorentz 
invariant. This gives a very different wave-functional from that given by the semi-classical 
expansion which would contain propagators, and so be non-local, but then we should not 
expect the strong-coupling limit to resemble the semi-classical one. To leading order the 
ground-state wave-functional is 

* [A] =exp-^ [ d 3 xtrB 2 , 7 2 = n 9*<W) N = ^ (1 1Q) 

Now that we have an approximation for \& we can compute expectation values of 
operators 0[A, E]. We will be principally interested in the Wilson loop, in which case O 
can be taken independent of E, so O = a; [A], say. (More generally we use the Schrodinger 
representation for E, (1.10) to reduce the action of n on $ to multiplication by some 
functional dependent on \I>.) 



_ /PANAMA] tt[A] 
{) ~ J VA **[A] \P[A] ' [ > 



and to leading order this is 



( DAw[A] exp ( A [ d 3 xtrB 2 ) 

<"> " —r TT^r V 1 ' (1.21) 

fVAexp(^fd 3 xtrB 2 ) 

But this is the functional integral for the vacuum expectation value of uj in three dimen- 
sional Yang-Mills theory, Wick rotated to Euclidean space, and with coupling 7. We are 
back to the problem we started with but in one dimension lower. Feynman has discussed 
the qualitative structure of the wave-functional for three-dimensional Yang-Mills theory, 
arguing that there is a mass gap [5]. We want to calculate (1.21) for large g, this means 
that we can repeat the above analysis. We begin by undoing the Wick rotation and working 
in three-dimensional Minkowski space with time coordinate x 3 , and quantisation surface 
x 3 = constant. We will denote vectors in the 12-plane by a bar. Working in the Weyl 
gauge A3 = the wave-functional of the vacuum is a functional of the two components of 
the gauge potential Ai,A2, and in the Schrodinger representation 



The Schrodinger equation is 



B = e ij (ViAj+AiAj) (1.23) 
as well as a constraint that demands that $ be gauge-invariant. 
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We again construct a Hamiltonian operator that has a finite action on a space of sums of 
products of local integrals with finite coefficents. Consequently the regulated Laplacian 
depends on a new arbitrary function, 0, and with B = —\ f d 2 xtrB 2 , Schrodinger's 
equation becomes 

[~\l 2 h + JS*) = (1-25) 
Expanding in a power series in g~ 8 



oo 



= exp _ = exp - £ ^7" 8 ^ n [A] , £ = 9 4 - 8n £n (1-26) 

\ n=l / n=l 

we obtain to leading order 

kA^S 1 =2(£ 1 -B/k). (1.27) 

Again acts on integrals of local gauge-invariant functions to generate further integrals 
of local gauge-invariant functions, and in the next section we compute 



A^B = fix (J d z xj +/3 2 B (1.28) 
With /3i = (N 2 - l)<p)(0)/87r, (3 2 = -23N/6n for SU(N). Taking 

2 * a h ( 



Si = -z-;B, Si = -f- / d'x) (1.29) 
P2K 2 fan \J J 

yields a solution to (1.27). The leading order expression for $[A] is thus 

*[A] = exp^ J d 2 xtrB 2 , ^ = (1.30) 

If we restrict our attention to those tv[A] that can be expressed as functionals of only A 
on the quantisation surface i.e. a; [A] = u[A] then 

so to leading order in the strong-coupling expansion 

^ fVA**[A]Q[A]*[A] ^ I vA ^ A ] exp{^Jd 2 xtrB^ 
fVA$*[A]$[A] JVAexp[^Jd 2 xtrB 2 ) ' 
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This is a functional integral for Yang-Mills theory in two Euclidean dimensions so we are 
back to where we started, but in two dimensions lower. We will now stop this process of 
reducing the dimension, because the gauge-invariant two-dimensional theory is free in the 
gauge A 2 = and we can evaluate (1.32) directly. 

It is well known that Yang-Mills theory confines in two dimensions in the sense that 
the Wilson loop has an area law behaviour. The Wilson loop is the vacuum expectation 
value of the trace of the path-ordered exponential of the gauge potential A taken around 
a closed curve [1]. If the curve is a rectangle of sides R and T then in the Euclidean 
formulation of the theory this has the interpretation of being exp — E(R)T, where E is 
the energy of a pair of massive 'test-quarks' created for a time T and held a distance R 
apart. With an area-law behaviour E(R)T = oRT so the potential energy of the two static 
quarks is oR, requiring an infinite energy to completely separate them. We will now show 
how the Wilson loop for four dimensional Yang-Mills theory shows area-law behaviour in 
the strong-coupling limit. In the Euclidean formulation of the theory we want to compute 

/ VA„ (tr P exp (- § c dx ■ A)) exp (-X J d 4 x tr F^F^) 

W[C] = — ^ r '- (1.33) 

/ VA„ exp (jij / d 4 x tr F^F^j 

P denotes path-ordering of the Lie algebra generators around the closed curve C. If C 
is planar we can use rotational invariance in four-space to place it in the 12-plane. Then 
W[C] is given by (1.20) with u = tr P exp — § c eix • A and to leading order by (1.32) with 
uj = tr P exp — § c dx • A. Evaluating this in the gauge A 2 = gives 

/ VAi (tr P exp (- § c dx 1 A x )) exp J d 2 x tr (d 2 A 1 ) 2 ) 



/ VA 1 exp / d 2 x tr (d 2 A l )^ 



tr P<.rp[-?± (j) (j) dx 1 dy 1 g(x,y)T A T A ) , (1.3 i) 



'C Jc 

where Q is the Green's function for <9|, i.e. 

d\Q{x,y) =8(x,y). (1.35) 
Bose symmetry requires that Q be symmetric in its arguments, so we take 

g(x,y) = ^\x 2 -y 2 \5(x 1 -y 1 ). (1.36) 

For the Lie algebra su(N) we have T A T A = —1(N 2 — 1)/N. Because this is pro- 
portional to the identity and because of the 5(x 1 — y 1 ) in the Green's function, the path- 
ordering is redundant. Furthermore the area enclosed by C can be written as 

A[C] = -<t I dx 1 dy x ^\x 2 - y 2 \8(x 1 - y 1 ). (1.37) 
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Putting all this together gives the leading order approximation for the Wilson loop as 



W[C] ~ Nexp (-*A[C\) , a = [ — J7 * = XN 2 (N 2 - 1)(0'(O))V . (1-38) 

Where A is a numerical constant, approximately equal to 2.6 x 10 -4 . Now a is a physical 
quantity which must not depend on the arbitrary function 0, so the coupling must actually 
be a functional of (p. <j> is dimensionless, but its argument has the dimension of [length], 
so 4>'(0) = \x has the dimension of [mass]. Since our expansion is in increasing powers of 
l/g 4 we shall define a (3- function 



(L39) 



where the differentiation is taken by holding physical quantities such as a fixed. Applying 
this to the above expression for a gives 

0=1 (1.40) 

which tells us that l/g 4 (/j,) decreases, as the mass scale is reduced (which we can also see 
directly from (1.38)). We should emphasise at this point that our expansion in powers of 
l/g 4 is of the ground-state wave-functional, and not the physical quantity a. Thus if we 
were to compute the first correction to (1.38) and use this to calculate the Wilson loop it 
would give corrections to a of higher order in l/g 4 , as well as, possibly, a contribution of 
the same order in g as our leading approximation, (1.38). It is inevitable in any quantum 
theory with only a single dimensionless coupling, where dimensionful physical quantities 
acquire those dimensions through dimensional transmutation that those quantities cannot 
be obtained from perturbation series in powers of the coupling, or powers of any function 
of the coupling, f(g) say. All physical masses, for example, must depend on / and the 
arbitrary mass scale in the same way, i.e. as [i exp — df / Pf(f'), with (3f = in order 
to be independent of \x. This would mean that if we could expand two physical masses 
in power series in / we would only need to compute the first terms in order to know the 
ratio of the two masses to all orders. Put another way if / — > as — > /U* say, then the 
/(/-independent ratio could be computed in this limit from only the first terms in the series. 
Such a situation could not take into account the full non-linearity of the theory. That our 
'strong-coupling' expansion of the ground-state wave- functional gives a non-zero result for 
a tells us that necessarily this will not be the first term in a power series expansion in the 
coupling. We will not be despondent about this, as the alternative is an approximation 
scheme in which we cannot compute a at all. However this may make the renormalisation 
of the expansion difficult. 

The wave-functional (1.19) is Lorentz invariant. Under a Lorentz boost the gauge con- 
dition Aq = must be preserved, so the change in A under an infinitesimal transformation 
has two pieces, the first being the usual one appropriate to the spatial components of a 
four-vector, whilst the second is a gauge transformation that preserves the Weyl gauge. 



9 



The generator is thus L(a) = ia • J d 3 xx.H, with H the Hamiltonian density. This can be 
regulated in the same way as H in the next section. Now, by construction, H^> = £^> , so 
Ti^f = (V -3 + const.)^, for some J. If we take / Ax = 0, then L(a)^> = -(/ d 3 ia-J)$. 
J has dimensions of [length] -3 , and to the order we are working this can only come from 
TC acting on B/fi, so J must be n~ x multiplied by a vector of dimension [length] -4 , but we 
cannot construct a local gauge-invariant vector with this dimension out of A and positive 
powers of derivatives of 0, so J must vanish, as can be seen by explicit calculation, so 
L(a)tf = 0. 

2. Heat-Kernel Regulator 

We will now address the crucial question of how to regulate Schrodinger's equation 
(1.11). The problem is that it contains a product of two functional derivatives at the 
same spatial point. So, for example A acting on — tr A 2 (x) is proportional to 5(0) which 
is meaningless. More generally, if E denotes the space of sums of products of three- 
dimensional integrals of local functions of A with finite coefficents, then A acting on a 
typical element of S does not yield another element of E. We will replace A in Schrodinger's 
equation by a regulated operator, A^, which does have a finite action on elements of 
E and maps E to itself. We will insist that for those a G E for which Act G E the 
action of A and A^ coincide. In addition we require that A^ is invariant under gauge 
transformations, as well as more general transformations of coordinates on configuration 
space, consistent with the interpretation of A as a Laplacian. This regularisation and 
renormalisation procedure introduces an arbitrary function on which A^ depends. This 
is an inevitable consequence of isolating the finite part of a divergent quantity, and has the 
effect of introducing an arbitrary scale into the problem. Ensuring that physical quantities 
are independent of this arbitrariness, e.g. by absorbing it all into the coupling, is the 
basic problem of renormalisation and has the effect of making the coupling a functional of 
4>. To construct A^ we begin by 'point-splitting' the functional derivatives. Consider the 
operator 

t ^/^^' kHS < x -^^ (2 - 1) 

Where the kernel, K, is a second rank tensor in three-space, satisfying some generalised 
heat equation ((2.16) below) so that it depends smoothly on the proper time t, with the 
initial condition 

Limit tiO K RS (x,y;t) = 8(x-y)5 RS l (2.2) 

it has the interpretation of being the 'temperature' at the point x at time t due to a 
delta-function source placed at y at the initial time. 1 is the identity tensor (dyadic) in 
three-space. Taking t small but non-zero, t = e say, gives a regulated operator that acts 
on elements of E to produce other elements of E. T e acting on a G E produces local 
integrals because in the small time e heat cannot flow vary far from its original source and 
so K(x, y, e) will depend only on A(x) for x in a neighbourhood of y of volume of order 
y/e . However, the coefficents of the integrals resulting from the action of T e on a G E will 
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contain powers of e some of which diverge as e — > 0, so we cannot simply replace Act by 
the Limit e — > T e cr. These powers of e may be determined from dimensional analysis. 
We will see that the proper-time e has dimensions of [length] 2 , whilst T has dimensions 
of [length] -1 . We will choose K to preserve gauge and rotational invariance, so, given the 
local character of T e , we have, for example 



/ 



d 3 x (aie 



T e B = 

" 5/2 + a 2 e~ 1/2 *rB 2 + a^'Hr ((VAB + AAB + BA A) 2 ) + ...) , (2.3) 



where the ai are numerical constants. We can extract a finite result for T e a when e = 
using analytic continuation in an analogous fashion to zeta- function regularisation [6]. 
Suppose 4>{p) is a differentiable function equal to one at the origin, and that f(p) is also 
a differentiable function which is finite at the origin and such that 4>(p)f(p) vanishes at 
infinity. Then on integrating by parts 



1(a) = / dpp fl -V(p)/(p) 
Jo 

For s > the first term vanishes, so 



-mm 

s 



dp P -^{H)- (2.4) 
s dp 



r°° d 

Limit sjO sl(s) = - dp — Uf) = f(0). (2.5) 

Jo d P 

If / diverges at the origin like p~ n , with n an integer, then /(0) is meaningless, but we can 
give a meaning to the left-hand side of (2.5) and use this as our definition of /(0). The 
integral I(s) exists (provided 0/ has no other divergences) for s > n, so if we analytically 
continue sl(s) to small values of s we can take the limit as s I to obtain a finite result. 
Thus, repeated integration by parts for s > n gives 

= D-d -^ WP>/WP " > ' <2 ' 6) 

This expression allows an analytic continuation to s < n so we can take the limit to obtain 

1 d n 

Limit siO sl(s) = -_(0/p")| p=0 . (2.7) 

nl dp n 

Clearly this result depends on the arbitrary function 0, as we should expect. Extracting 
a finite quantity from a divergent one requires the use of an arbitrary renormalisation 
procedure in general and the introduction of an arbitrary scale in particular. Physical 
quantities will be independent of (p if the theory is renormalisable. With the use of this 
tool we take the action of on an arbitrary element, a G E to be 



/>oo 

A^cr = Limit s I si dpp s ~ 1 (fi(p)T p 2cr, 

Jo 



(2.8) 



where the integral on the right-hand side is defined by analytic continuation in s as in the 
example above. We take T at proper-time p 2 rather than p because of the square-roots 
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appearing in (2.3) which are generic to a theory in an odd number of dimensions, (by 
which we mean here the three spatial dimensions). For example 



A^B = y>z(^0( 5 )(O) + a 2 0'(O)£rB 2 ). (2.9) 

More generally, acting on a local integral of dimension [length] ~ D results in other local 
integrals of dimensions [length] ~ D with D' < D. We have done more than regulate A 
by constructing an operator whose action on elements of S does not depend on a cut-off. 
We have actually done part of the job of renormalisation. The rest, ensuring that physical 
quantities are independent of <p, remains to be done. 

We now turn to the choice of heat-kernel. Under an infinitesimal gauge transformation 
parametrised by u> 

A c (x) -> A c (x) + V^ c (x) - f RSC co R (x)A s (x.), (2.10) 
the functional derivative transforms homogeneously as 

6 6 /^V)ttIt^ (2-11) 



5A c (x) 5A c (x) J v J 5A s (x) 

so that A is invariant. We will demand that A^ is also gauge-invariant, which in turn 
requires that T t be gauge-invariant, so we insist that the kernel transform under (2.10) as 



K^(x, y; t) -> K flS (x, y; t) -f AB V(x)K BS (x, y; t) -K* B (x, y; t)f ABS u A (y). (2.12) 
Or in matrix notation with uj = (u RS ) = (u® f® RS ), 

K(x, y; t) -> K(x, y; t) + w(x) K(x, y; t) - K(x, y; t) to(y). (2.13) 

It is easy to construct a kernel that satisfies these conditions. If A is a local operator that 
acts on variations of the gauge-potential as 

(Ao5A) c (x) = J d 3 y A CD (^y)-dA D (y), (2.14) 

then we can view A as a matrix with indices x, y, CD as well as the three-space tensor 
indices, and o as matrix multiplication. If, in addition, A transforms under (2.10) as 

A(x, y) - A(x, y) + w (x) A(x, y) - A(x, y) u(y), (2.15) 

then we can take the kernel to be exp — tAo acting on the identity matrix. For example, if 
V is the covariant derivative acting in the adjoint representation , i.e. VidAj = VidAj + 
[A t , 5Aj] then we could (but we won't) take A o 5 A = —ViViSA, in which case A(x, y) = 
T>iDi5(x., y) satisfies (2.15). We can then compute K by solving the heat equation 

d 

—A o K = — K, (2.16) 

dt v ' 
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subject to the initial condition (2.2). 



There are many choices of operator A that satisfy the above conditions. We need a 
further criterion to pick the right one. This will emerge from considering the geometry that 
underlies Schrodinger's equation (1.11). If we think of A(x) as being a set of coordinates 
on the space of configurations, A, then this space has a natural (flat) metric, 

g{8A,8A) = - J d 3 xtr8A(x) -<5A(x), (2.17) 

which is invariant under (2.10). We can use this metric to define the volume element T>A 
in (1.20). Furthermore the operation o is just the scalar product constructed using this 
metric. The operator A is then the Laplacian, expressed in terms of the coordinates A, 
acting on functions on A such as i.e. if D/DA is the Levi-Civita connection on A 
(indistinguishable from 8/8 A for the metric (2.11)) then 

i / D D \ D D . . 

A = 9 {da'da) = da°da' < 2 ' 18 > 

and 

D oKo-^-. (2.19) 



DA DA 

Physics should be independent of how we choose coordinates on A, so Schrodinger's equa- 
tion must transform covariantly under a change of coordinates A(x). \I/ transforms as 
a scalar, and A is the scalar Laplacian so B = — \ J d 3 xtr~B 2 must also transform as a 
scalar. Now we want the regulated operator T to act on scalars to produce scalars, and 
this will be the case if K is a second rank tensor under transformations of A. This in turn 
will be guaranteed if A is also a second rank tensor. So we need to find a local operator in 
three-space enabling us to construct the heat equation (2.10), which is also a tensor under 
coordinate transformations of A. We can construct tensors by differentiating scalars using 
the Levi-Civita connection. The only scalar we have identified so far is £>, so 

D 8 

B= — B = VAB + AAB + BA A = VAB (2.20) 



DA 8A 

is a covariant vector on A, whilst a suitable candidate for A is 

With the flat metric (2.11) this gives 

A o 8 A = -V ■ V8A + VV ■ 8 A + 2B A 8 A + 28 A A B 



= —V ■ V8A + W ■ 8 A — 2F • 8 A (2.22) 

This is just the operator that appears in the semi-classical expansion of three-dimensional 
Yang-Mills. F is the field-strength tensor in the adjoint representation. 
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We need to compute the coefficents a±, a 2 that occur in the expansion (16). Now 

tb= Ja ( oKo (is 8 )) = {jK oK ) o JA. B + Tr { Ko JK {jK B )) ' (2,23) 

Tr denotes a trace over all the matrix indices. The first term on the right-hand side is 
proportional to DAB and so contributes to 013 but not ol\ and «2- We are left with 

Tr ( K ja (tz b ) ) = Tr (K A) = ~it Tr (K) (2 - 24) 

Evaluated at t = e. So we need to compute the diagonal elements of K. We could attempt 
to do this for a general A, but since we only need to extract the numerical coefficents 
«i, «2 we can do the calculation for a background satisfying the classical equation of 
motion D A B = 0. We will show that in this case 

exp (-tAo) = exp (t(V ■ VI + 2F)-) - I dr V exp (rV ■ V) V- (2.25) 

Jo 

(K can be found for a general configuration by perturbing about this result in powers of 
I>AB.) For an arbitrary A 

[V-V,V] = -2F-V + (V-F), (2.26) 

so that 

AoVuj = -{V -F)uj, D-Ao 5 A = (V -F) ■ 5 A. (2.27) 

From now on we will assume that A satisfies T> A B = which implies that T> ■ F vanishes 
so then Vu is a zero-mode of Ao and V ■ (A o 6 A) = 0. Let V denote the projector onto 
the kernel of V- so that 

Vo5A = 5A-V(V-V)~ 1 V-5A, (2.28) 

Now A o V = A so 

exp (-tAo) VoSA = exp (-tAo) 5A - V (V ■ D)" 1 V ■ 5 A, (2.29) 
If u is in the kernel of V- then 

Aou = -D-Du-2F-u (2.30) 
which will also be in the kernel of V- so we can write 

exp (-tAo) V = exp (t(V ■ VI + 2F)-) V (2.31) 

From (2.26) 
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(V-V1 + 2F) -V = VV-V, (2.32) 

so 



exp (t(V ■ VI + 2F)-) V = Vexp (W ■ V) , (2.33) 

(as can be seen by expanding both exponentials in power series, and using (2.32) to move 
V to the left.) This enables us to simplify (2.31) to 



exp (-tAo)Vo = exp (t(V -VI + 2F)-) - Vexp (W ■ V) (V ■ V) 1 V- (2.34) 
Finally we can write (2.29) as 

exp (-tAo) =exp (t(V • VI + 2F)-) — Vexp (W-V) (V ■ V)' 1 V ■ +V (V-V)~ 1 V- 

(2.35) 

which is equivalent to (2.25). Using cyclicity of Tr we obtain 

—TrK= — (Trexp (t(V ■ VI + 2F) ■) — Tr exp (W-V)). (2.36) 
We can make this more explicit by introducing the heat-kernels 

H CD (x, y, t) = exp (W ■ V) CD 5(x - y), 

G CD (x, y, t) ee exp (t(V ■ VI + 2F)-) CD 6(x - y) 1, (2.37) 
which satisfy the generalised heat equations 

d d 

—H = VWH, — G = (V ■ VI + 2F) • G, (2.38) 
subject to the initial conditions 

Limit t[0 H BC (x,y) = 5(x-y)5 BC , 

Limit tiO G sc (x,y) =5(x-y)5 BC l. (2.39) 
In terms of these functions the right-hand side of (2.36) becomes 



d 3 x G£ c (x, x, t) - / d s x # cc (x, x, t) , (2.40) 



d_ 
dt 

so we need to know the diagonal matrix elements of G and H for short times. These may 
be computed using a standard technique [7] that assumes the existence of expansions of 
the form 
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_( x _y) 2 /4i 00 

_( X _y) 2 /4( OO 

G BC (x, y, f) = 6 -=3 £ ^ b n c (x, y). (2.41) 

V47Tt n=0 

The prefactor solves the ordinary heat equation and collapses to a delta-function at small 
times, consequently the initial conditions are satisfied by taking 

a^ c (x, x) = 5 BC , b* c (x, x) = 5 BC 1. (2.42) 

Substituting the expansions into the heat-equations (2.38) and equating terms of the same 
order in t leads to recurrence relations 

(x-y)-Db (x,y) = 
(V 2 1 + 2F) • b n _! - (x - y) • T?b n (x, y) = nb n (x, y), n > 
(x-y) -£>a (x,y) =0, 

V 2 a n -i - (x - y) • X>a„(x,y) = na n (x,y), n > 0. (2.43) 

Repeated applications of V to these equations, followed by setting x = y enables them to 
be solved for b n (x, x), a n (x, x) and their derivatives. For example 

ai(x,x) = 0, a 2 (x,x) = —FijFij, 

bi (x, x) = 2F, b 2 (x, x) = 2F • F + -L 1 F l3 F l3 + ^V-VF (2.44) 

Using these in (2.40) and (2.23) gives eventually 

3n„ 11C , , 

ax = ^3, a 2 = 3 (2-45) 

V4n 6\/47T 

where n g is the number of Lie algebra generators and C is given by f Rp Qf s Q p = —§ RS (J. 
For SU(N) we have n g = N 2 — 1 and C = 2N. Note that if we had performed this 
calculation for an Abelian theory, such as QED, then A would not depend on A and 
consequently the right-hand side of (2.9) would just be independent of A. We would then 
be unable to solve (1.16) with our locality ansatz. 

To construct the Schrodinger equation for three-dimensional Yang-Mills theory we 
need an operator that acts on an element, a, of a space E of two-dimensional integrals 
of local functions of A to produce other elements of that space. Proceeding as before 
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we define an operator A^ so that its action on a coincides with that of the unregulated 
Laplacian when the latter is finite. Thus 

/•OO 

A 7(7 = Limit s I s dpp s ~ 1 (p(p)T p a, 

Jo 

To ensure gauge and reparametrisation invariance the kernel is chosen as K = exp — tAo 
with A a tensor in the space of gauge-potentials, 

A RS (x, y) = - - g. (2.47) 

As before, K can be expressed in terms of simpler heat-kernels 

K(x, y) = G(x, y) - I drVH(x, y)V-, (2.48) 
Jo 

with Vf = Vf+ [A, f] and 

—H = V-VH, —G=(V-m + 2F)-G, (2.49) 

Here Fij = eijB a dj is the field strength in the adjoint representation. G and H are assumed 
to have small t expansions of the form 

p -(x-y) 2 /4t °° 

G BC (x^t) = ^ t X>"ff g (x,y), 

n=0 
— (x— y) 2 /4t 00 

tf BC (z, y, f ) = 6 4?rt £ raf c (x, y). (2.50) 

n=0 



The calculation of the expansion coefficents is the same as in the previous case, leading to 

1 

6" 



ai(x,x) = 0, a 2 (x,x) = -B 2 adj , 



11 „0 =r 1 

"6 

Using these expansions we find 



h(x, x) = 2F, b 2 (x, x) = ~ — B 2 adj l + -V-VF (2.51) 



T e B = J d 2 x (&! e" 2 + a 2 tr B 2 + ...) (2.52) 



with a\ = a 2 = |ff , from which follows (1.28). 
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3. Conclusions 



We have constructed a regulated form of Schrodinger's equation for four-dimensional 
continuum Yang-Mills theory in which the Hamiltonian operator has a finite action on 
integrals of local functions of A. We solved this for the logarithm of the ground-state 
wave- functional as a power series in l/g 4 on the assumption that the terms in this series 
are themselves integrals of local functions of A. Ground-state expectation values calculated 
using the leading order result reduce to calculations in two-dimensional Yang-Mills theory. 
This is known to confine and be directly related to certain string theories. Consequently 
the Wilson loop in the four-dimensional theory satisfies an area law, and the coupling is 
seen to grow as the mass-scale, n, is reduced, which is infra-red slavery. It is the intrinsic 
non-linearity of Yang-Mills theory that is responsible for this. We would not have been able 
to solve Schrodinger's equation for an Abelian gauge theory using the ansatz of locality. 
Our results are preliminary, and many questions remain to be answered. We need to know 
how to ensure that all physical quantities are independent of the arbitrary function that 
was introduced to regulate, and partially renormalise the theory. We need to compute 
corrections, and to construct the excited states. We also need to know how to incorporate 
fermions. We expect that other quantum theories such as CP n -models and String Theories 
may be treated in a similar fashion. 

Finally, it is a pleasure to acknowledge conversations with Ed Corrigan on the subjects of 
zeta-function regularisation and Yang-Mills theory, and also with Ryu Sasaki. 
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